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ABSTRACT 

The brightest supernova remnant in the Magellanic Clouds, N132D, belongs to the rare class of 
oxygen-rich remnants, about a dozen objects that show optical emission from pure heavy-element 
ejecta. They originate in explosions of massive stars that produce large amounts of O, although 
only a tiny fraction of that O is found to emit at optical wavelengths. We report the detection 
of substantial amounts of O at X-ray wavelengths in a recent 100 ks Chandra ACIS observation of 
N132D. A comparison between subarcsecond-resolution Chandra and Hubble images reveals a good 
match between clumpy X-ray and optically emitting ejecta on large (but not small) scales. Ejecta 
spectra are dominated by strong lines of He- and H-like O; they exhibit substantial spatial variations 
partially caused by patchy absorption within the LMC. Because optical ejecta are concentrated in a 
5 pc radius elliptical expanding shell, the detected ejecta X-ray emission also originates in this shell. 
Subject headings: ISM: individual (N132D) — supernova remnants — X-rays: ISM — supernovae: 
general 



1. INTRODUCTION 

The heavy-element ejecta of core-collapse (CC) super- 
novae (SNe) are dominated by oxygen. By studying the 
properties of O and other freshly synthesized elements 
in supernova remnants (SNRs), we can infer the total 
amount of ejected O, relative abundances, and the spa- 
tial distribution of heavy elements. With this informa- 
tion, we can then constrain the progenitor main-sequence 
mass and learn more about how massive stars explode. 
Oxygen is most readily seen in X-rays, as evidenced by 
a growing number of O-rich SNRs detected by modern 
X-ray satellites such as Chandra and XMM-Newton. X- 
ray detection of plasma enriched in O but deficient in Fc 
allows us to identify CC SNRs in the first place. High- 
resolution studies of SNRs containing both optical and 
X-ray-emitting O ejecta hold particular promise for ad- 
vancing our understanding of ejecta in CC SNRs. 

N132D, an X-ray-bright SNR in the Large Magel- 
lanic Cloud ( LMC), contains optic a lly emitting O-rich 



ejecta ( Sutherland fc Dopita||1995| |Morse et al-111995. 
19961. Morse et al. (1995) interpreted their observations 



m terms of an expanding (1650 km s _1 ) 10 pc ellip 
soidal ejecta shell. Assuming undecelerated expansion, 
they estimated the remnant's age at 3150 y r 10 times 
ol der tha n Cas A. Spectrosco py with ASCA (|Hughes et 
aT] |1998[ ) and XMM-Newton ( |Behar et aLpOOip showed 
that the X-ray emission is dominated by the shocked am- 
bient medium. X-ray emission from ejecta has not yet 
been convincingly demonstrated, although it might be 
expected based on the high (30-35 M ) mass of the SN 
progenitor inferred from optical and UV spectroscopy 
(Blair et al. 20001. We report here the discovery of 
clumpy X-ray-emitting O ejecta in N132D with new 
Chandra X-ray observations, matching (on large scales) 
the optical ejecta morphology as seen by Hubble. 
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N132D was observed by the Chandra X-Ray Observa- 
tory on 2006 January 9, 10, and 16 with the Advanced 
CCD Imaging Spectrometer (ACIS) S3 CCD chip, at the 
same telescope roll angle and target location for each ob- 
servation, for a total effective exposure of 89.3 ks. Data 
were processed with CIAO version 3.4 and CALDB 3.3.0 
with default processing options. While the data were ac- 
quired in the very faint telemetry mode, we chose the 
default faint format in data analysis because of a signif- 
icant loss of source photons in high surface brightness 
sections of the remnant with the very faint format. A 
large background region covering most of the S3 chip 
was used for all spectra. Spe ctral analysis was per- 
formed with XSPEC version 12 ( IArnaud|1996| . We used 
the nonequilibrium ionization (NEI) version 2.0 ther- 
mal models, based on the APEC/APED spectral codes 
( Smith et al.| | 200TJ) and augme nted by addition of inner- 
shell processes (Badcncs et al. 2006). 

We obtained 2.8 x 10° source photons, with about 0.5% 
background. Images were smoothed using platelets with 
a default smoothing parameter 7 = 1/2 (Willett 2007); 
Figure [l] (top) shows the remnant's X-ray morphology. 
We measured the position of a bright star visible in X- 
rays on the S3 chip to be within O'.'l of its USNO UCAC2 
position. 

High spatial resolution imaging of O-rich ejecta was 
done by the Hubble Space Telescope W ide Field Planetary 
Camera 2 (WFPC2; [Morse et al.| 19961 and the A dvanced 
Camera for Surveys (ACS; Beasley et al. |2004[ ). Hubble 
ACS images obtained on 2004 January 22 are most use- 
ful, because they provide a complete spatial coverage of 
the remnant, and the broad bandpass of the ACS Wide 
Field Channel (WFC) F475W filter includes [O III] A5007 
emission from O ejecta even with most extreme radial 
velocities. In addition to the ACS F475W data, we also 
use ACS images taken with the F658N and F775W fil- 
ters, to highlight the shocked ambient medium emitting 
strongly in Ha and to distinguish diffuse SNR emission 
from numerous LMC stars. Individual ACS exposures 
(four per filter) were aligned and combined, and cosmic 
rays were removed with the on-the-fly reprocessing sys- 
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Chandra ACIS image o f N132D in th e 0.3-7.0 keV 
energy range, smoothed with platelets (Willctt 2007). The image 
is 120" X 115" in size. North is up and east is to the left. Note 
the complex filamentary structure within the remnant's interior 
(enclosed by an ellipse) and the bright outer blast wave. Bottom: 
A three-color Hubble ACS image in F475W, F658N, and F775W 
filters (in red, green, and blue, respectively). O-rich ejecta are 
prominent in the F475W filter, while the shocked ISM radiates 
predominantly in the F658N filter (Beasley et al. 2004; see also a 
high quality image in Brown 2007, p. 26). The location of optically 
emitting O-rich ejecta is marked in both images. 

tern at STScI. We refined the astrometric accuracy of the 
final combined images by identifying and measuring posi- 
tions for a number of Guide Star Catalog II (GSC2) stars 
and shifting images from the GSC1 reference frame to 
the more accurate (0'.'3) GSC2 frame. We present these 
images as a three-color image in Figure [I] ( bottom) ; red 
diffuse emission marks the location of O-rich ejecta. We 
outline optically emitting O-rich eje cta, detected both in 



Hubble and in ground-based data of |Morse et al.] ( 1995[ ), 
with solid lines in the bottom and top panels of Figure 

A similar approach allows us to locate O-rich ejecta 
in X-rays. We constructed a composite three-color im- 
age from soft (0.3-0.5 keV), medium (0.5-0.75 keV), and 




Fig. 2. — Merged image between 0.3 and 7 keV. Red: 0.3- 
0.5 keV; green: 0.5-0.75 keV; b lue: 0-75-7 ke V. All three images 
were smoothed using platelets (Willctt 2007|. X-ray knots and 
filaments near optical O-rich ejecta are shown in the inset. Red 
ellipses show ejecta knots and filaments with particularly strong 
emission in the O band (0.5-0.75 keV). Regions in green contain a 
mixture of ejecta and the shocked ambient medium. Three selected 
blast wave regions are enclosed by red rectangles. 

hard (0.75-7 keV) energy bands (Fig. [2j) . Emission in the 
soft band is particularly sensitive to variations in the in- 
terstellar medium (ISM) absorption toward N132D. The 
medium band contains strong O Hea 0.57 keV and O 
Lya 0.65 keV lines; O-rich gas will produce an excess of 
emission in this band when compared with the ambient 
LMC gas. The broad hard band measures the strength of 
the overall X-ray emission. Figure [2] reveals widespread 
variations in X-ray colors in the central region. Enhance- 
ments in the medium band (greenish knots and filaments) 
correlate spatially with optical O-rich ejecta, implying 
O-rich X-ray-emitting ejecta with a similar large-scale 
distribution as the optical ejecta. 

We attempted to quantify this spectral distinction be- 
tween knots near optical O-rich ejecta and emission else- 
where in the interior. We selected a large number of X- 
ray emission enhancements in the N132D interior (within 
the ellipse shown in the top panel of Fig. [T] chosen to 
exclude the bright outer blastwave). This was done by 
repeated thresholding of the image shown in Figure [l] 
(top), filtered first by the unsharp masking procedure, 
and further processed with the morphology-based oper- 
ation "opening" (a basic image-processing algorithm) to 
spatially separate partly overlapping emission enhance- 
ments. This procedure provided us with many spatial re- 
gions, preferentially sampling X-ray knots and filaments, 
as opposed to more spatially uniform inter-knot emission. 
Our final sample includes 145 regions. 

We then counted photons within each region in un- 
smoothed Chandra images in each of the three energy 
bands, arriving at soft (S), medium (M), and hard 
(H) counts. We first examined regions not near O- 
rich ejecta, by excluding the regions shown in Figure [l] 
In Figure pi we plot their colors C$ = log(S/H) arid 
Cm = log(M/H). We also excluded regions with less 
than 16 counts (corresponding to a signal-to-noise ratio 
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Fig. 3.— Color-color scatter plot, 0.3-0.5/0.75-7 vs. 0.5- 
0.75/0.75-7 keV count ratios (on log-log scale), for central regions 
of N132D (within the ellipse of Fig. 1, top). Colors of X-ray- 
emitting gas not overlapping spatially with O-rich optical ejecta 
cluster along a linear regression line with slope 3/2, while red and 
green regions associated with optical ejecta (Fig. 2) exhibit emis- 
sion excess in the 0.5—0.75 keV band. Red circles denote regions 
with the most extreme colors, more than 3<r away from the regres- 
sion line. Three outer blast wave regions are marked by diamonds. 

S/N < 4) in the soft band, in order to reduce effects 
of counting (Poisson) noise. There is a large spread in 
X-ray colors in the center of N132D (Fig. [3]), but they 
are strongly correlated. A linear regression describes this 
correlation well, and scatter about the regression line is 
dominated by the counting statistics in the soft band. 

We next selected 31 X-ray emission enhancements lo- 
cated in close proximity to O-rich ejecta; they are shown 
in Figure [2j Their location in the color-color plot differs 
from the other central regions. While their colors are 
also strongly correlated, they show excess emission in the 
medium band (Fig. [3]) . Regions with the most extreme 
colors appear green in Figure [2j we marked them by large 
circles in Figure [3] The remaining regions appear blue; 
Figure[3]reveals that they still show an excess emission in 
the medium band. This excess emission indicates X-ray- 
emitting O-rich ejecta. The spatial distributions of op- 
tical and X-ray-emitting ejecta are strongly correlated, 
although a detailed comparison of Figures fl]( bottom) and 
[2] reveals that there is no exact correspondence between 
optical and X-ray emission on the smallest spatial scales. 
Clumpy X-ray-emitting ejecta have the same large-scale 
distribution as opticall y emitting ejecta, f orming an ex- 
panding elliptical shell ( Morse et al.|1995 1 . 



3. SPECTROSCOPY 

In order to investigate the large spread in X-ray colors 
in N132D, we examined spectra of the blast wave at a 
number of locations along its periphery. We found sub- 
stantial spatial variations in the ISM absorption, which 
greatly affect spectra at low photon energies. Three blast 
wave spectra in Figure [4] demonstrate th is effect. We 
mo deled them with a plane shock model (Borkowski et 
al. 2001 1, using a two-component (Galactic plus LMC) 
absorption. The Galactic absorption is 5.5 x 10 20 cm - 
at the N132D location (based on H I radio observations 
of Staveley-Smith et al.|2003[ ), while the intrinsic absorp- 
tion JMh within the LMC was allowed to vary. We as- 
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Fig. 4. — Spectra of the blastwave in three locations shown in 
Fig. 2, one in the west (crosses, red) and two in the south (stars, 
green; and filled circles, blue). Plane-shock model fits are shown 
by solid lines. Large spectral variations at low energies are caused 
by spatially varying absorption within the LMC. 

sumed 0.4 solar abundances within absorbing LMC ma- 
terial, except for nitrogen, whose abundance is low (0.1 
solar) in the LMC. The same abundances have been as- 
sumed for the emitting gas, except for Ne, Mg, and Fe 
(with Ni tied to Fe), whose abundances were allowed to 
vary. The best-fit value for N H (LMC) of 1.4 x 10 20 
cm~ 2 ([0, 4] 90% confidence interval) in the western blast 
wave region is consistent with no absorption within the 
LMC, while significant absorption, 1.6(1.2,2.1) x 10 21 
cm~ 2 and 4.1(3.6,4.7) x 10 21 cm~ 2 is required in the 
two southern regions. Plasma temperatures are nearly 
equal, 0.66, 0.71, and 0.67 keV, respectively — while ion- 
ization ages are longer in the south [6.1(3.8,8.4) x 10 11 
cm~ 3 s and 15(8.7, 23) x 10 11 cm' 3 s] than in the west 
[4.0(2.9,5.3) x 10 11 cm' 3 s]. 

The X-ray colors of the blastwave depend strongly on 
Nff] colors corresponding to spectra of Figure |4j and 
plotted in Figure [3] spread along the regression line. Col- 
ors become bluer down the line as absorption increases. 
We conclude that spatial variations in interstellar absorp- 
tion are mostly responsible for the large range in X-ray 
colors seen in Figure [2] The significant difference in Nh 
between the two adjacent regions in the south demon- 
strates that absorption is clumpy on small spatial scales, 
although a larg e-scale north - south absorption gradient 
is also present. Banas et al. (1997) argued that N132D 
is physically associated with a molecular cloud located 
~ south of the remnant. While there is litt le spatial 
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pact Array H I 21 cm survey images ( Kim et al. 2003 1 



reveal the presence of substantial amounts of H I emis- 
sion around this molecular cloud, extending to the north- 
ern tip of N132D. It is likely that this gas just outside 
the molecular cloud is responsible for the patchy X-ray 
absorption. 

Spectra of several regions associated with optically 
emitting O-rich ejecta and with excess emission in the 
medium band are shown in Figure [5] These regions are 
marked by thick red ellipses in Figure [2] we added spec- 
tra for adjacent regions if their colors were similar. The 
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Fig. 5. — X-ray spectra of ejecta in four locations shown in Fig. 2 
by thick red ellipses (we added spectra from adjacent ellipses in 
order to improve S/N ratio). For comparison, an X-ray spectrum of 
the shocked ambient gas is also shown (Fig. 2, thick green ellipse). 
From top to bottom (at 0.6 keV): North, east, ambient, west, and 
central spectra, with colors (Cm, Cs) equal to (0.15, -0.80), (-0.10, 
-1.14), (-0.64, -1.62), (-0.17, -1.34), and (-0.47, -1.72), respectively. 
Prominent X-ray lines are labeled. Oxygen lines in ejecta spectra 
are stronger than in the blast wave spectrum. 

top spectrum corresponds to two northern knots with ex- 
tremely soft hardness ratios (these knots are located in 
the top right corner of Figure [3| . Bright optic ally emit- 
tin g O-r ich ejecta filaments overlap spatially (Morse et 



1995[ ) at the location of the northern knots. Strong 
mes dominate the X-ray spectrum; they are much 



stronger than O lines in the blast wave spectra of Figure 
[4| The O Hea line of the He-like +6 ion is stronger 
than the O Lya line of the hydrogenic +7 ion. It is 
likely that a lower than average ionization state of O 
accounts for the softness of the X-ray spectrum. Spec- 
tra of ejecta in the east and west, corresponding to th e 



blueshifted filaments B4 and Bl of Morse et al. (1995), 



also show strong O lines. These are most typical of ejecta 
spectra in N132D. The O Lya line is now stronger than 
the O Hea line, resulting in harder spectra and bluer X- 
ray colors. The X-ray colors of the eastern knots place 
them in the main "clump" of ejecta knots in the X-ray 
color-color plot (Fig. [3]), significantly below the color- 
color relationship for the shocked ambient gas. The X- 
ray spectrum of a small central ejecta knot (the bottom 
spectrum in Fig. 51) corresponds to even bluer X-ray col- 
ors (Cm = —0.47, Cs — —1.72). A comparison with 
the hard (Cm = —0.64, Cs = —1.61) spectrum of a 
bright interior knot with the normal (LMC) abundance 
(Fig. [5| reveals the presence of stronger than average O 
lines, Dut they are weaker than in other ejecta spectra, 
likely because of the increased absorption. A substan- 
tial amount of absorbing LMC gas appears to be present 
toward the very center of N132D; a central "blue hole" 
seen in Figure [2] corresponds to an absorbing ISM cloud 
with A^(LMC) = 3 x 10 21 cm' 2 . 

Strong O lines in ejecta spectra may be due either 
to enhanced O abundances or to short ionization ages. 
We fit the top spectrum of Figure [5] with a plane shock 
model with temperature of 0.59 keV and ionization age 
of 4 x 10 10 cm -3 s and with nearly the same abundances 
as in the blast wave. For electron densities of > 15 cm -3 



(typical of the N132D blast wave; [Williams et~aT|[ 2006 ) , 
the very short ionization age implies a shock age of less 
than 100 yr. This is clearly an unreasonably short shock 
age for a 3000 yr old SNR. Problems arise also if low 
(~ 1 cm~ 3 ) electron densities are assumed, as such low- 
density gas with LMC abundances is too faint in X-rays. 
This would also require a substantial (an order of mag- 
nitude or more) pressure imbalance between the ejecta 
and the blast wave, rather difficult to sustain for a pro- 
longed period of time. These severe difficulties suggest 
that the alternative explanation of strong O lines in terms 
of O-rich ejecta is the correct one. The determination 
of plasma conditions and abundances within the ejecta 
would require a separate, detailed spectroscopic analy- 
sis involving a superposition of heavy-element ejecta and 
blast wave emission. We just note here that the low O 
Lya/O Hea line ratio of 0.41 (0.36 after correction for 
absorption) in the plane shock fit just mentioned implies 
a low (0.2 keV) plasma temperature under conditions 
of collisional ionization equilibrium expected in dense 
shocked O-rich ejecta clumps. It is likely the observed 
clumpy ejecta emission arises in relatively dense O-rich 
ejecta, which have been shocked by low-velocity shocks 
propagating into O-rich clumps. We estimate that a few 
tenths of M of O is present in these clumps. 

4. DISCUSSION 

High spatial resolution imaging and spectroscopy of 
N132D with Chandra reveal the presence of clumpy O- 
rich ejecta in its center. The spatial distribution of X- 
ray-emitting ejecta correlates well with optically emit- 
ting ejecta on large (but not small) scales. This shows 
that optical and X-ray emission can be closely linked in 
dynamically advanced remnants such as N132D, where 
most of the O might have already been shocked. In 
N132D, ejecta are located in an expanding ellipsoidal 
shell. No such shell is present in either G292.0+1.8 or 
Pup A, two O-rich Galactic SNRs with ages similar to 
N132D. Optically emitting ejecta in G292.0+1.8 are dis- 
tributed throughout much of th e remnant, with little de - 
tailed correlation with X-rays ( Winkler & Long 2006 1 . 
Asymmetrically distributed optical ejecta filaments in 
Pup A are moving in the oppo site direction to its neutron 



star (Winkler & Petre 2007). This demonstrates that 
the innermost metal-rich ejecta are strongly affected by 
poorly understood processes leading to CC explosions. 
No correlation of optical emission with more evenly dis- 
tributed X-ray-emitting ejecta has been reported. X- 
ray expansion veloc ities in N132D are poorly known; 
Hwang et al. (19931 found (marginal) evidence for an 
X-ray expansion consistent with the optical expansion 
in spatially integrated spectra obtained by Einstein, but 
archival Chandra and XMM-Newton grating spectra are 
more suitable for studying the ejecta kinematics. The 
origin of the X-ray and optically emitting ejecta shell is 
unknown at this time. It may be the location of a reverse 
shock propagating into ejecta; more likely, it may be an 
ejecta shell created by a Ni bubble effect shortly after the 
SN explosio n, similar to what is s een in the O-rich SNR 
B0049-73.6 flHendrick et~aT||2005j ) . More detailed anal- 
ysis of existing Chandra and X MM-. 
is warranted. 
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